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Abstract 



breifly described. It makes possible to analytically investigate and calculate 



m 

\ A new zero modes enhancement (ZME) model of the true QCD vacuum is 

O ' 
CO 

I numerically low-energy QCD structure from first principles. Expressions of 

in ' 

0^ ' basic chiral QCD parameters (the pion decay constant, the quark and gluon 



uum energy density (up to the sign, by definition, the bag constant), suitable 



Qh| condensates, the dynamically generated quark mass, etc) as well as the vac 

P 

for numerical calculations, have been derived. Solution to the Schwinger- 
Dyson (SD) equation for the quark propagator in the infrared (IR) domain 
on the basis of the ZME effect in QCD was used for this purpose. There are 
only two independent quantities (free parameters) by means of which calcu- 
lations must be done within our approach. The first one is the integration 
constant of the above mentioned quark SD equation of motion. The second 

one is a scale at which nonperturbative effects begin to play a dominant role. 
PACS numbers: 11.30 Rd, 12.38.-t, 12.38 Lg and 13.20 Cz. 
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I. THE ZERO MODES ENHACEMENT MODEL OF QUARK CONFINEMENT 

AND DCSB 



Today there are no doubts left that the dynamical mechanisms of quark confinement [1] 
and dynamical chiral symmetry breaking (DCSB) [2, 3] are closely related to the complicated 
topological structure of the QCD nonperturbative vacuum. For this reason any correct 
nonperturbative model of quark confinement and DCSB necessary becomes a model of a 
true QCD vacuum and vice versa. Also, it becomes clear that nonperturbative infrared (IR) 
divergenses are closely related to the above mentioned nontrivial vacuum structure, play an 
important role in the large scale behaviour of QCD [4-6]. If it is true that QCD is an IR 
unstable theory (i.e., if has no IR stable fixed point) then the low-frequency modes of the 
Yang-Mills fields should be enhanced due to nonperturbative IR divergences. So the full 
gluon propagator can diverge faster than the free one at small momentum, in accordance 
with [7-10] 

d,M) ~ (q 2 r 2 , q 2 - (1.1) 

which describes the zero modes enhancement (ZME) effect in QCD. If, indeed the low- 
frequency components of the virtual fields in the true vacuum have a larger amplitude than 
those of the bare (perturbative) vacuum [5], then the Green's function for a single quark 
should be reconstructed on the basis of this effect. It is important to understand that a 
possible effect of the ZME (1.1) is our primary dynamical assumption. We will consider this 
effect as a rather realistic confining ansatz for the full gluon propagator in order to use it as 
input information for the quark SD equation. 

Let us consider the exact, unrenormalized SD equation for the quark propagator in 
momentum space (Fig. 1) 

S-\ V ) = S \p) + g 2 C F J ^^(p, q)S(p - qhvD^(q), (1.2) 

where Cf is the eigenvalue of the quadratic Casimir operator in the fundamental represen- 
tation. Since the notions are obvious, let us only note that our parametrization of the full 
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quark propagator is as follows 

- t S(p)=pA(-p 2 ) + B(-p 2 ). (1.3) 
The full gluon propagator in the arbitrary covariant gauge is 
D^(q) = -i 



q 2 



l A - q \a) + a q -^ , (1.4) 



q' q H 



where a is a gauge fixing parameter (a = for Landau gauge). 
Assuming that in the IR region 

d(-q 2 , a) = + (3(a) + 0(q 2 ), q 2 -> 0, (1.5) 

where n is the appropriate mass scale parameter, we obtain the above mentioned generally 
accepted form of the IR singular asymptotics for the full gluon propagator (1.1) [4, 5, 7-10] 
(enhancement of the zero modes). 

In order to actually define an initial SD equation (1.2) in the IR region (at small mo- 
menta) let us apply the gauge-invariant dimensional regularization method of 't Hooft and 
Veltman [11] in the limit n = 4 + 2e, where e is a small IR regulation parameter, which is to 
be set to zero at the end of the computation e — > + . Because of this the quark propagator 
and other Green's functions become dependent, in general, on this IR regulation parameter 
e. We consider the SD equations and the corresponding quark-gluon ST identity in Eu- 
clidean space (d n q — > id n q E , q 2 — > —q%, p 2 — ► ~Pe, but for simplicity the Euclidean 
subscript will be omitted). 

Let us use, in the sense of the distribution theory, the relation [12] 

(q 2 y 2+€ = ^ 4 (q) + {q 2 )' + 2 + 0(e), e - + , (1.6) 

where (q 2 )+ 2 is a functional acting on the main (test) functions according to the so-called 
"plus prescription" standard formulae [12]. Substituting (1.4-1.6) into the quark SD equa- 
tion (1.2) on account of the above mentioned "plus prescription" formulae and expanding 
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in powers of q and keeping the terms of order q~ 2 ( the Coulomb order terms), in agree- 
ment with (1.5), one finally arrives at the quark propagator expansion in the IR region for 
unrenormalized quantities [13] (in four-dimensional Euclidean space). 

As mentioned above, all Green's functions become dependent generally on the IR reg- 
ularization parameter e. In order to extract the finite Green's functions in the IR region, 
we introduce the renormalized (IR finite) quark-gluon vertex function at zero momentum 
transfer and the quark propagator as follows T M (p, 0) = Zi(e)r /X (p, 0) and S(p) = Z 2 (e)S(p) 
at e ^ + , respectively. Here Z^e) (i = 1,2) are the corresponding IR renormalization 
constants. The e-parameter dependence is indicated explicitly to distinguish them from 
the usual UV renormalization constants. f 0) and S(p) are the renormalized (IR finite) 
Green's functions and therefore do not depend on e in the e — > + limit, i.e. they exist as 
e — > + . The correct treatment of such strong singularity (1.6) within distribution theory 
[12] enabled us to extract the required class of test functions in the renormalized quark 
SD equation. The test functions do consist of the quark propagator and the corresponding 
quark-gluon vertex function. By the renormalization program we have found the regular (at 
zero point) solutions for the quark propagator (see below). For that very reason relation 
(1.6) is justified, it is multiplied by the appropriate smooth test functions. Due to a quark 
convergence condition, a cancellation of nonperturbative IR divergences takes place. Be- 
cause of this condition the explicitly gauge-dependent terms (the so-called next-to-leading 
terms) in the above mentioned quark propagator expansion become e - order terms. For this 
reason these noninvariant terms vanish in the e — > + limit [13]. 

Absolutely in the same way should be reconstructed the ghost self-energy and the corre- 
sponding ST identity for the quark-gluon vertex [13]. We develop a method for the extraction 
of the IR-finite Green's functions in QCD which means that they do not depend on the IR 
regulation parameter e as e — > + . For this purpose we have worked out a renormalization 
program in order to cancel all the IR nonperturbative divergences which makes it possible to 
obtain a close set of the SD equations and the corresponding ST identity in the quark sector. 
By completing our renormalization program, we explicitly show that only multiplication by 
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the quark IR renormalization constant will remove all nonperturbative IR divergences from 
the theory on a general ground. We have shown [13] that for the covariant gauges the 
complications due to ghost contributions can be considerable in our approach. 

The closed set of equations, obtained in our paper [13], which will be used for numerical 
calculation of basic chiral QCD parameters as well as the vacuum energy density, should 
read 



where Sq 1 (p) is the free quark propagator. The finite at zero point the ghost self-energy is 
denoted by 6(0) and g 2 includes the mass scale parameter /j 2 , determining the validity of the 
above mentioned deep IR singular asymptotic behaviour of the full gluon propagator (1.5). 
Let us mention that the IR finite quark renormalization constant, explicitly not shown here, 
which should multiply the free quark propagator in Eq. (1.7) is to be set to unity without 
losing generality (multiplicative renormalizability) [13]. It is worth noting also that Eq. 
(1.7) and Eq. (1.8) describe the leading terms of the corresponding expansions of the quark 
SD equation and ST identity in the IR region, respectively [13]. 

In order to solve the system (1.7-1.8), it is convenient to represent the quark-gluon vertex 
function at zero momentum transfer as follows 



S- 1 (p) = S 1 (p) + ~g 2 T,(p,0)S(p) ll 



(1.7) 



^6(o)f>, o) = id,s-\ P ) - ±b(o)s(p)r,( P , o)s-\p) 



(1.8) 



r„(p, 0) = F^p 2 )^ + F 2 (p 2 ) PfM + F 3 (p 2 )p,p + F 4 (p 2 )p7, 



(1.9) 



Substituting this representation into the ST identity (1.8), one obtains 



-A( P 2 ), 



F2{p 2 ) 



-2B\p 2 ) - F,(p 2 ) 



F 3 (P 2 ) 



F,(p 2 ) 



2A\p 2 ), 
A 2 {p 2 )B-\p 2 ) 
E{p 2 ) 



(1.10) 
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Here the prime denotes differentiation with respect to the Euclidean momentum variable p 2 
and A(p 2 ) = A(p 2 )E-\p 2 ), B(p 2 ) = B(p 2 )E-\p 2 ) with E(p 2 ) = p 2 A 2 (p 2 ) + B 2 (p 2 ). For 
the sake of convenience, the ghost self-energy at zero point b = 6(0) is included into the 
definition of a new coupling constant A = g 2 \b(0)]' 1 (2n)^ 2 . 

Proceeding now to the dimensionless variables by p 2 = n 2 t = [i 2 \z and parameters 
jt = z ,t = and introducing then the following notations 

A(p 2 ) = ^ 2 A{t) = -^ 2 jg(z), B 2 {p 2 ) = ^ 2 B 2 {t) = fi- 2 jB 2 (z , z), (1.11) 

and doing some algebra, the initial system (1.7-1.8) can be rewritten as follows (normal 
form) 

g'(z) = -[- + l]g(z) + - + m B(z) (1.12) 

z z 

B'(z) = ~g\z)B-\z) - [m g(z) + B(z)}, (1.13) 

where rho = mo(f) 1 / 2 . It is easy to check that solutions to this system in the chiral limit 
nio = are: 

g( z ) = z' 2 [exp(-z) - 1 + z] (1.14) 

and 

zo 

B 2 (z , z) = 3 exp (-2z) J exp (2z')g 2 (z') dz . (1.15) 

2 

The exact solutions (1.14) for g(z) and (1.15) for the dynamically generated quark mass 
function B(z ,z) are not entire functions. The functions g(z) and B(z ,z) have removable 
singularities at zero. In addition, the dynamically generated quark mass function B(zo,z) 
also has algebraic branch points at z = z and at infinity. Apparently, these unphysical 
singularities are due to ghost contributions. The quark propagator may or may not be an 
entire function but in either cases the pole-type singularities should disappear. This is a 
general feature of quark confinement and holds in any gauge. 
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In order to reproduce automatically the correct behaviour of the dynamically generated 
quark mass function at infinity, it is necessary to put z = oo in (1.15) from the very 
beginning. Obviously, in this case solution (1.15) cannot be accepted at zero z = 0, so one 
needs to keep the constant of integration z arbitrary but finite in order to obtain a regular, 
finite solution at zero point. 

The region z > z can be considered as nonperturbative, whereas the region z < z can 
be considered as perturbative. By approximating the full gluon propagator by its deep IR 
asymptotics such as (q 2 )~ 2 in the whole range [0, oo), we nevertheless obtain that our solution 
for the dynamical quark mass function B{zq,z) manifests the existence of the boundary 
value momentum (dimensionless) z which separates the IR (nonperturbative) region from 
the intermediate and UV (perturbative) regions. If QCD confines then a characteristic scale, 
at which confinement and other nonperturbative effects become essential, must exist. On the 
other hand, because of this one can eliminate the influence of the above mentioned unphysical 
singularities, coming from the solutions to the quark SD equations (due to necessary ghost 
contributions), on the S'-matrix elements reproducing physical quantities. 

Thus, within our approach to QCD at large distances [13] in order to obtain numerical 
values of any physical quantity, e.g. the pion decay constant (see below), the integration over 
the whole range [0, oo] reduces to the integration over the nonperturbative region [0, z ] , 
which determines the range of validity of the deep IR asymptotics (1.1) of the full gluon 
propagator and consequently the range of validity of the corresponding solutions (1.14) and 
(1.15) for the IR piece of the full quark propagator. 

Let us make the main conclusions now. First, if the enhancement of the zero modes of the 
vacuum fluctuations (1.1) takes place indeed then the quark Green's function, reconstructed 
on the basis of this effect, has no poles. In other words, the enhancement of the zero modes 
at the expense of the virtual gluons alone removes single quark from the mass-shell (quark 
confinement theorem of Ref. 13). Second, a chiral symmetry violating part of the quark 
propagator in this case is automatically generated. From the obtained system (1.12-1.13) 
it explicitly follows that a chiral symmetry preserving solution (m = B{z) = 0,g(z) ^ 0) 
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is forbidden. So a chiral symmetry violating solution (m = 0,B(z) ^ 0,g(z) ^ 0) for the 
quark SD equation is required. Thus the enhancement of the zero modes automatically 
leads to quark confinement and DCSB at the fundamental quark level and they are in close 
connection with each other. 

II. THE VACUUM ENERGY DENSITY AND GLUON CONDENSATE 

The effective potential method for composite operators [14] enables us to investigate 
the vacuum of QCD since, in the absence of the external sources, the effective potential is 
nothing but the vacuum energy density, the main characteristic of the vacuum. The effective 
potential at one-loop level is [14] 



where S(p) (1.3), S (p) and D{p) (1.4), D (p) are the full, free quark and gluon propagators, 
respectively. Here and everywhere below the trace over space-time and color group indices 
is assumed but they are suppressed. Let us recall that the free gluon propagator can be 
obtained from (1.4) by setting simply d(—q 2 , a) = 1. 

Evidently the effective potential is normalized as follows ^(S'o, -Do) = V(Sq) — V(D ) = 
0. Because of this normalization the vacuum energy density now should be defined as 
follows e = V(S ,D ) — V(S,D) = —V(S,D). This means that the perturbative vacuum is 
normalized to zero and e = e q + e g with e q = —V(S), e g = —V(D), where V(S) and V(D) 
are given by (2.1). 

Going over to Euclidean space (d A p — * id n p, p 2 — > — p 2 ) and dimensionless variables 
and parameters (1.11), we finally obtain after some algebra (n = 4) 



V(S, D) = V(S) + V(D) 




(2.1) 



3 




/ 



dzz{lnz zg 2 (z) + B 2 (z ,z) - 2zg(z) + 2}, 



(2.2) 



8tt 2 



o 
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where we introduced the UV cutoff which should be identified with the arbitrary constant 
of integration z as was discussed in previous section. The explicit expressions for g(z) and 
B 2 (zo,z) are given by (1.14) and (1.15) respectively. The expression (2.2) describes the 
comtributions to the vacuum energy density of the light confining quarks wuth dynamically 
generated masses. 

The vacuum energy density due to the nonperturbative gluon contributions in the same 
variables is 



€ 9 



= lA; V / dz z {ln(l + *-) - A + b} (2.3) 
o 

Here one important remark is in order. In fact, vacuum energy density e g does not vanish 
at z — > oo as it should because of the above mentioned normalization. Thus it needs an 
additional regularization at this limit. From (2.3), it follows that the term containing the 
constant b should be substracted from this expression. So regularized vacuum energy density 
should be calculated through the relation (2.3) which becomes 



f - Li- 4 *- 2 x 

t 9 ~ ff2 



181n(l + |)-i, 2 ln(l + ^)-^ ' 



(2.4) 



The vacuum energy density due to confining quarks (2.2) automatically disappears at z — > 
oo, so it does not requires any additional regularization. 

The vacuum energy density is important on its own right as the main characteristics of 
the nonperturbative vacuum of QCD. On the other hand, it makes possible to estimate such 
important phenomenogical parameter as the gluon condensate introduced within the QCD 
sum rules approach to resonance physics [15]. Indeed, through the vacuum energy density e 
it can be expressed as follows 

(Ol^GyO) = -f c = -f (e g + e 9 ). (2.5) 

The weakness of this derivation [15] is, of course, that it was obtained on the basis of 
the perturbative calculation of the f3{a s ) -function. In any case, it would be instructive to 
estimate the gluon condensate with the help of (2.5). 
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III. BASIC CHIRAL QCD PARAMETERS 



Because of its especially small mass, the pion is the most striking example of the Gold- 
stone realization of chiral symmetry SU (2) L x SU (2)r. Beside the quark condensate and the 
dynamically generated quark mass it is one of the three important chiral QCD parameters 
that determine the scale of chiral dynamics. 

I. The pion decay constant F w is defined in the current algebra (CA) as 

(014(0)1^(9)) =<F^. (3.1) 

(The normalization F n = 92.42 MeV is used [16]). Clearly, this matrix element can be 
written in terms of the pion- quark-antiquark proper vertex and quark propagators as 

iF^F = J (04^{(|) W*S(p + q)Qi(p + q,p)S(p)}. (3.2) 

To get expression for F n one has to differentiate Eq. (3.2) with respect to the external 
momentum q v and then set q = (Fig. 2). 

Information on the BS pion wave function up to terms of order q can be obtained from 
the corresponding axial-vector vertex. Indeed, in our paper [17] it has been found that this 
vertex can be decomposed in a self-consistent way into pole (dynamical) and regular parts 
as follows 

r l 5 > + q,p) = -^Gl(p + q,p) + rg(p + q,p), (3.3) 
where the BS bound-state amplitude is 

Gl(p + q,p) = ~ (jj ^G(p+q,p), (3.4) 

with 

G(p + q,p) = 2B(-p 2 ) + qR 6 (-p 2 ) +pqR u (- p 2 ) (3.5) 

and arbitrary form factors are the residues of the corresponding form factors, entering the 
vertex from the very beginning [17]. The regular part at zero momentum transfer q = is 
determined as 
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rg(p,p) = (yhsb^i +p^G 2 +p^G 3 +p^G 4 } 



(3.6) 



where in the Euclidean space 



-A(p 2 ) - i? 6 (p 2 ) 

2EV) 

2A'(p 2 ) 



G 2 (p 2 ) 



G 3 (P 2 ) 



G 4 (p 2 ) 



i?n(p 2 ). 



(3.7) 



This system is nothing else but the conditions for the cancellation of the dynamical singu- 
larities at q = for the corresponding form factors [17]. The regular part at q = also 
depends on the same form factors i? 6 and Ru as the BS bound-state amplitude up to terms 
of order q. 

Then by taking into account the BS pion wave function up to terms of order q, as given 
by (3.4-3.5), and (1.3) with the substitution p — > p + q and expanding in powers of q, the 
expression (3.2) can easily be evaluated. Going over to Euclidean space (d 4 p — > id 4 p, p 2 — > 
— p 2 ) and using dimensional variables (1.11), one finally obtains 



Here the primes denote differentiation with respect to the dimensionless Euclidean momen- 
tum variable t, A — A(t), B = B(t) and quantities with overline are shown after Eqs. (1.10). 
Here and in what follows we denote the pion decay constant F n in the chiral limit of the CA 
representation by Fqa- 

The main problem now is to find a good nonperturbative ansatz for both arbitrary form 
factors Rj(t)(J = 6, 11) in the IR region, which determine the contributions to the pion decay 
constant from the second diagram in Fig. 2. In the nonperturbative calculations these terms 
cannot be ignored by saying formally they are of order g 2 in the coupling constant, as it 
was done in the perturbative treatments [3, 18-19]. In connection with this let us point 




(3.8) 
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out that the difference between the vector and axial-vector currents disappears in the chiral 
limit. For this reason let us assume [20] that the IR finite quark-gluon vector vertex function 
at zero momentum transfer (1.9-1.10) is a good approximation to the regular piece of the 
axial- vector vertex at zero momentum transfer in the chiral limit (3.6-3.7). A fortunate 
feature that admits to exploit partial analogy between vector and axial-vector currents in 
the chiral limit for the flavor non-singlet channel is that the contribution to the pion decay 
constant in the CA representation (3.8) does not depend on the form factor G 2 (p 2 ) at all. 
In this case the analogy between (1.10) and (3.7) becomes complete, and one obtains 

W)=0, Rn(p 2 ) = -^^^. (3.9) 

Of course, we cannot prove these relations but it will be shown later (Part II) that this 
dynamical assumption (nonperturbative ansatz) leads to very good numerical results for all 
chiral QCD parameters thereby justifying it once more. 

In terms of the new parameters and variables (1.11) and on account of (3.9) we finally 
recast (3.8) as follows 

F CA = ^kW J q dz ( 3.io) 

This expression will be used for numerical calculation of the pion decay constant in the CA 
representation. 

II. The other important chiral QCD parameter is the dynamically generated quark mass 
rrid, defined as the inverse of the full quark propagator (1.3) in the chiral limit at zero point 
[13, 17, 21] 

m d =[iS ch (0)}-\ (3.11) 

where S c h(0) denotes the full quark propagator in the chiral limit m = 0. Obviously, this 
definition assumes also regularity at the zero point. Though the dynamical quark mass rrid 
is not an experimentally observable quantity, by all means, it is desirable to find such kind 
of solutions to the quark SD equations in which dependence on a gauge-fixing parameter 
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disappears. In this sense rrid defined by (3.11) becomes gauge- invariant. As it was briefly 
discussed in Section 1, exactly such a nonperturbative quark propagator has been found 
within our approach to QCD at large distances [13]. 

Using the standard decomposition of the quark propagator (1.3) and its inverse, dynam- 
ical chiral symmetry breaking (DCSB) at the fundamental (microscopic) quark level can be 
implemented by the following condition 

{S- 1 (p), l5 } + = il,2B(- P 2 ) ± 0, (3.12) 

so that the 75 invariance of the quark propagator is broken and the measure of this breakdown 
is the double of the dynamically generated quark mass function 2B(—p 2 ). Precisely this 
quantity at zero 2B(0) can be defined as a scale of DCSB at the fundamental quark level 
[17]. In accordance with (3.11), let us denote it by 

A CSBq = 25(0) = 2m d , (3.13) 

The definitions rrid and A C sBq have now direct physical sense within the above mentioned 
solutions to the quark SD equation. 

Let us write down the final result for the dynamically generated nonperturbative quark 
mass (3.11) too, expressed in terms of the new parameters and variables (1.11) 

m d = k {z B 2 (z ,0)}- 1/2 , (3.14) 

where B 2 (zo,0) is given by (1.15) at zero point. 

III. As it is well known, the order parameter of DCSB - quark condensate can also be 
expressed in terms of the quark propagator scalar function B(—p 2 ) (1.3). Its definition is 

(qq) = (0|gg|0> = / ^JrS(p). (3.15) 

The final result expressed in terms of new variables and parameters (1.11) is as follows 

_ 3 _3/2 f z ° 

{QQ)o = ~^2 k o z o J dzzB(z ,z), (3.16) 
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wher the superscript denotes the chiral limit. As a function of rrid the quark condensate 
can be expressed on account of (3.14). 

Thus there are only two independent (free) quantities by means of which all calculations 
should be done in our approach. The first one is the constant of integration z of the 
above mentioned quark SD equation of motion. The second quantity is a scale k at which 
nonperturbative effects begin to play a dominant role. The ZME model of a true QCD 
vacuum enables us to describe quark confinement and DCSB on a general ground. We begin 
Part II with numerical investigation of the low energy QCD structure at the chiral limit. 
At low energies QCD is under control of SU^Nf) x SUr{Nj) chiral symmetry (Nf is the 
number of different flavors) and its dynamical breakdown in the vacuum to the corresponding 
vectorial subgroup. So to correctly calculate basic low energy QCD parameters in the chiral 
limit means to corectly understand the dynamical structure of QCD at low energies. That 
is why it is important to start from the chiral limit. 
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FIGURES 



FIG. 1. The SD equation for the quark propagator in momentum space. 
FIG. 2. An exact expression for the pion decay constant F n in the current algebra (CA) 
representation. Here G J 5 , S and are the pion-quark bound-state wave function, the quark 
propagator and the axial-vector current, respectively. The slash denotes differentiation with respect 
to momentum q u and setting q = 0. 
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